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Abstract We derive the neutrino flavor transition probabilities with the neutrino treated as a wave packet. The 
decoherence and dispersion effects from the wave-packet treatment show up as damping and phase-shifting of the 
plane-wave neutrino oscillation patterns. If the energy uncertainty in the initial neutrino wave packet is larger than 
around 0.01 of the neutrino energy, the decoherence and dispersion effects would degrade the sensitivity of reactor 
neutrino experiments to mass hierarchy measurement to lower than 3 a confidence level. 
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1 Introduction 

Much information regarding neutrino mixing have been revealed in the past few decades. However, in most oscillation 
data analyses, neutrinos are described as plane waves with definite energy and momentum. Since neutrino produc¬ 
tion and detection are spatially localized, a wave-packet description is more general and appropriate for a complete 
understanding of neutrino oscillations. Even if the plane-wave treatment is a good approximation for neutrino flavor 
transitions, the wave-packet decoherence and dispersion effects could still give rise to small corrections to oscillation 
parameters. We investigate the wave-packet treatment in detail, constrain the energy uncertainty er W p of reactor anti¬ 
neutrinos, and calculate corrections to the mixing parameters by the Daya Bay [Tj and KamLAND [ 2 ] reactor neutrino 
experiments. 

According to our analyses, the wave-packet treatment does not produce significant modifications of the mixing 
parameters measured by current reactor neutrino experiments based on the plane-wave analysis. However, current 
experimental data allows a large possible range in the initial momentum width of the neutrino wave packet (ov). 
If the initial momentum/energy uncertainty of the neutrino wave packet is larger than around 0.02 of the neutrino 
energy, the decoherence and dispersion effects could have significant effects on future measurements of the neutrino 
mass hierarchy. 

In this article, we apply a wave-packet treatment to neutrino oscillations El H El H 0 El 0 and examine its 
phenomenological implications on reactor neutrino experiments at medium baseline. 

The article is organized as follows: In Section 2, the wave-packet treatment is briefly reviewed. In Section 3, we 
use the survival probability derived in Section 2 to explore the sensitivity of potential measurements of mass hierarchy 
with medium baseline reactor experiments. We summarize and conclude in Section 4. 


2 Wave-packet impact on current reactor neutrino experiments 

2.1 Wave-packet treatment for neutrino oscillations 

The plane-wave description of neutrino oscillation has been developed for almost 40 years m- In the standard cal¬ 
culation of plane-wave neutrino oscillations, the probability of detecting a neutrino flavor state \vp) with energy E , 
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Fig. 1: An illustration of the decoherence and dispersion effects on neutrino flavor oscillations. The different mass 
eigenstates (red and blue lines) in a flavor neutrino would gradually separate in space due to speed differences, leading 
to decoherence and reduction of interference. At the same time, each mass eigenstate wave packet would disperse, 
partially compensating for the decoherence effect. However, the overlapping fractions of the wave packets remain 
smaller than those of plane waves. 


evolved from a pure flavor state \v a ), at a distance z from the production point is given by 

P a p = S a p - 4 ^2 Re(U a jUpjU* k Up k )sm 2 (Arnljz/4:E) 

k>j 

+ 2 ^ lm(U aj Up j U* k Up k )sm(Aml j z/2E), (1) 

k>j 

where U a i denote the elements of the PMNS matrix and Am 2 j = — m 2 are the differences of the mass eigenvalues 

squared mi- 

However, as neutrino production and detection are spatially localized, there must be finite intrinsic energy/momentum 
uncertainties and a neutrino should be described by a wave packet. The wave-packet character of light has been 
discussed in details in Ref. m- Based on similar arguments as in that reference, all particles are produced and de¬ 
tected as wave packets. In 1976, the wave-packet nature of propagating neutrinos was proposed m ■ A wave-packet 
description is expected to be more general and appropriate for a complete understanding of neutrino oscillations 
Nevertheless, there are also arguments against the wave-packet treatment. Refs. i TSl fl9] argue 
that a wave-packet description is unnecessary as the oscillation system is stationary. However, it has been pointed out 
in Refs. [201-21,, "221 [231 [23] that the authors of Refs. [181119] have mixed the macroscopic stationarity with microscopic 
stationarity. The wave-packet description of neutrino oscillations is necessary at least in principle. 

Therefore, a neutrino is described by a wave packet as it propagates freely mmm- 

(p-p»r 

2 erg 

( 2 ) 
( 3 ) 


kOM)) = f 


dp 


exp 


sj ^/iro, 
■exp [i(pz - Ei(p)t)]\vi), 
\v a {z,t)) = ^2u*i\vi(z,t)), 


where \vi) is an energy eigenstate with energy Ei, p v is the mean momentum, ay is the width of the wave packet in 
momentum space 1 , assumed to be independent of the neutrino energy here, and \v a ) is a neutrino flavor state. Fig. [l] 
pictorially describes the wave-packet effects on the propagations of neutrino mass eigenstates. 

In order to calculate the integral in Eq. we expand the energy Ei(jp) around the mean momentum p u up to 
second order 

TT\r o 

Ei(jp) « Ei(p„) + Vi(p v ){p - p„) + 2(Ei(p u )) 3 ^ P ~ P ^ * ( 4 ) 


1 Here, a u is the effective uncertainty, with l/cr^, = l/<7p rod + l/^ det , which has included both the production and detection neutrino 
energy uncertainties 11611211 [20) . Moreover, we would like to point out that cr det represents the energy uncertainty of detection at the 
microscopic level, ie., that of the inverse-beta decay reaction. This is different from the detector energy resolution, which is determined 
by macroscopic parameters such as the performance of PMTs and geometry of the anti-neutrino detector, etc. In principle, the detector 
resolution is irrelevant for the size of the neutrino wave packets. 
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where Vi(p v ) 
the neutrino 


dEi 

dp 


= p v IEi(jp v \ is the group velocity of wave packet. We use Eq. (2), (3) and (4) to calculate 


flavor transition probabilities at baseline L\ 
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Detailed derivation of Eq. © is shown in Appendix A. In Eq. the terms in the first bracket correspond to the 
standard plane-wave oscillation probabilities, and those in the second bracket represent the modifications from the 
wave-packet treatment. The exp(— Xij) term corresponds to the decoherence effect due to the fact that different mass 
states propagate at different speeds Vi(p„) and they gradually separate and stop to interfere with each other, resulting 
in a damping of oscillations. The terms depending on yij come from the quadratic correction in Eq. 0 ; they describe 
the dispersion effects and are dependent on the dispersion length(s) 2 Lfj s . Furthermore, yij are proportional to (Twp> 
while Xij oc <r wp only. Therefore, if <r W p •C 1, the dispersion effect is expected to be more suppressed and negligible. 
Dispersion has two effects on the oscillations. On the one hand the spreading of the wave packet compensates for the 
spatial separation of the mass states, hence restoring parts of their interferences. On the other hand, dispersion reduces 
the overlapping fraction of the wave packets, and thus the interference or oscillation effects cannot be fully restored. 
Moreover, it also modifies the flavor oscillation phases: 


, 2t tL (1 „i, , , \ 

faj = £osc + \^2^ an _ J 


( 6 ) 


with deviations from the standard plane-wave oscillation phase written in the parentheses. If yij = 0, then fcj just 
reduce to the standard plane-wave oscillation phases. 

Additional discussions about the phenomenological consequences of the wave-packet treatment and details of the 
derivation of Eq. © can be found in [25l |3 0 [201 EO • 

In this paper, we focus on the analyses of reactor neturino experiments. According to our wave-packet treatment, 
in reactor neutrino experiments, the anti-electron neutrino survival probability is 


Pee = 1 — ^COS 4 (013 )sin 2 (2012 ) [1 
isin 2 ( 20 i 3 )cos 2 ( 0 i 2 )[l 
isin 2 ( 20 i 3 )sin 2 ( 0 i 2 )[l 


1 + vh 

( _ 1 _ 
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)*exp(—A2i)cos(</>2i)] — 
) I exp(—A 3 1 )cos(0 3 i )] - 
) i exp(—A 3 2 )cos(0 3 2)]. 


( 7 ) 


2.2 The constraints from current reactor experiments 

To date, the value of the parameter <T wp has not yet been determined. If o- wp is not negligible, the wave-packet effects 
could be significant and have important implications on current and future neutrino oscillation experiments. In this 
article, we constrain er wp by analyzing the published Daya Bay and KamLAND data shown in references [26] and [2], 
considering statistical errors only. Figs. [ 2 ] and [ 3 ] show the data points from Daya Bay experiment [26], along with the 
oscillation curves corresponding to different values of a wp . Fig. [4] shows the data points from KamLAND [2] and the 
oscillation curves of plane-wave and wave-packet treatments. In Figs. [2] to [4] L e g- ect ive are the flux-weighted average 
reactor baselines 3 , and the error bars just show the statistical uncertainties. 

2 The “dispersion length” in this report represents the distance where the dispersion effect becomes important in the neutrino oscillation. 
A different definition of dispersion length and more detailed discussion of dispersion can be found in ESMZHXZ1. 

3 In Daya Bay, the effective baselines are calculated for all three experimental halls. 
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We use Eq. 0 to fit the data points in Figs. [ 2 ] and [ 3 ] to get the constraint of cr wp from the Daya Bay Experiment. 
In these figures, the black (solid) curves correspond to <r wp = 0, sin 2 20i3 = 0.084 and Ato § 2 = 2.43 x 10 - ' 5 eV 2 , same 
as those of the standard plane-wave treatment. The brown (dashed) curve represents the wave-packet result with <r wp 
= 0.1, sin 2 20i3 = 0.084 and Am§ 2 = 2.4 x 10 -3 eV 2 (the best-fit mixing parameters when er wp = 0.1). The green 
(dashed) curve shows the oscillation pattern with cr wp = 0.3, sin 2 2$i3 = 0.096 and Am ^ 2 = 2.22 x 10~ 3 eV 2 (the 
best-fit mixing parameters when cr wp = 0.3). The blue (dot-dashed) curve corresponds to cr wp = 0.5, sin 2 2$i3 = 0.108 
and Ato§ 2 = 2.10 x 10~ 3 eV 2 (the best-fit mixing parameters when <r wp = 0.5). The data points in Figs. [ 2 ] and [ 3 ] show 
clearly the existence of neutrino oscillation, and the data agree with the P ss derived by the plane-wave approach (Eq. 
0 ) for a certain set of mixing parameters. This is not surprising since the baselines of the Daya Bay Experiment are 
short compared to the coherence length, so that the oscillations are not washed out yet. 

Moreover, the decoherence and dispersion effects are dependent on the baseline L, and thus the wave-packet impact 
in the far hall of Daya Bay is expected to be more significant than in the near halls. In the plots in Fig. [3j the black 
(solid) and brown (dashed) curves overlap almost completely, while the differences between the black, green and blue 
curves are also small compared to the error bars. However, for the far hall EH3, the wave-packet impact becomes 
significant for cr wp > 0.3. It is because in Eq. 0 . the term A ij depends on the baseline L. Since the effective baseline 
of EH3 is longer, the damping of oscillation (decoherence effect) in EFI3 is more significant. 

The result of our data analysis is shown in Fig. [5} Here, we have just considered the statistical errors. The constraints 
on the parameters could become worse with the systematic errors taken into account. Fig.[5]shows that the wave-packet 
impact is not significant in Daya Bay experiment and it hardly affects the measurement of $ 13 . The vertical black 
line in the figure represents the best-fit value of sin 2 2 $i 3 in the plane-wave model. A larger value of er wp , or larger 
decoherence effect, implies that the true value of $13 should be larger 4 . From our analysis, there is no strong evidence 
to suggest non-zero er wp . Moreover, our result suggests that the modification of @13 is not significant even when the 
wave-packet framework is considered. The best-fit value of sin 2 2 ( 9 i 3 from plane-wave analysis (vertical black line) is 
not ruled out even with 1 a C.L. We believe that it is because the effective baselines of the Daya Bay Experiment are 
short compared to the coherence length. 

We perform the similar analysis with KamLAND data from Fig. [4j The result of our data analysis is shown in Fig. 
[ 6 ] Again, the systematic errors are not considered. 



Fig. 2: Pee as a function of L e g ec ti ve /E for Daya Bay EH3, compared with data |26| . The black solid curve is the 
standard plane-wave result (sin 2 20i3 = 0.084, Am§ 2 = 2.43 x 10 -3 eV 2 , cr wp is negligible), while the brown, green 
and blue curves are wave-packet results with cr wp = 0.1 (and the corresponding best-fit values of sin 2 2$i3 = 0.084, 
Ato§ 2 = 2.4 x 10“ 3 eV 2 ), 0.3 (and the corresponding best-fit values of sin 2 20i3 = 0.96, Ato§ 2 = 2.22 x 10 -3 eV 2 ) and 
0.5 (and the corresponding best-fit values of sin 2 2(9i3 = 0.108, Am 2 2 = 2.10 x 10 -3 eV 2 ) respectively, ^effective was 
obtained by equating the actual flux to an effective antineutrino flux using a single baseline. 


Eq. ([ 5 ]) shows that L coh /L dls = 2er wp , which implies that the dispersion length is much longer than coherence 
length if cr wp is negligible. This means that the decoherence effect from separation of wave packets is expected to be 
more significant than the dispersion effect. However, the dispersion effect is also important because it partly restores 

4 If <T wp is non-negligible but we still see oscillation effect from the data, it means that the true value of the mixing angle is actually 
larger than we expected in plane-wave assumption. 
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Fig. 3: Same as Fig. [ 2 ] but for Daya Bay EH1 (upper) and EH2 (lower). 



Fig. 4: Pee as a function of L e g ect i ve /E in the KamLAND experiment. The data points are the ratios of the geo¬ 
neutrino and background subtracted i7 e spectrum to the expected no-oscillation spectrum [2j. -^effective = 180 km is 
the flux-weighted average reactor baseline. The blue (solid) curve is the best-fit plane-wave result, with sin 2 2$i2 = 
0.857 and Am 2 1 = 7.53 x 10 -5 eV 2 ; while the red (dashed) curve is the best-fit wave-packet result, with er wp = 0.12, 
sin 2 2$i2 = 0.902 and = 8.06 x 10 -6 eV 2 . The error bars show the statistical uncertainties only. 


the oscillation, especially in the case of large <7 wp . If the dispersion effect is not considered, the modifications of the 
true values of mixing angles in Figs. [5] and [6] would be more significant. The bounds in Figs. [5] and [6] come from the 
combination of decoherence and dispersion effects, but the contribution from decoherence is expected to be dominant. 

Figs. [5] and [6] show that wave-packet effects are not significant in the current reactor neutrino experiments. The 
1 a upper bound of the energy uncertainty <r W p ~ 0(0.1), which is larger than some previous theoretical estimations 
((7 wp ~ O(10 -7 )) [27i l6l l8l |28j . Although our analyses have not considered systematic errors, our study on the current 
reactor experiments suggest that er wp can be around 0(0.1) for reactor experiments. We emphasize that a W p in this 
article is an effective parameter which include both the production and detection neutrino energy uncertainties. The 
estimation of the value of this parameter or the size of neutrino wave packet has not come to a strong conclusion 
yet. Conventionally, it has been argued that a wp should be much smaller than 1 [28]• However, there are still no 
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Fig. 5: The 1 <j (solid), 2 a (dashed) and 3 a (dotdashed) constraints on “er wp vs sin 2 20i3” from fitting the published 
Daya Bay data [26]. Only statistical errors are considered. We have marginalized over Am |i in this plot. The vertical 
black line represents the best-fit sin 2 20 i 3 based on the plane-wave analysis. 



sin 2 20i2 

Fig. 6: The 1 a (solid), 2 a (dashed) and 3 er (dotdashed) allowed regions on “<r W p vs sin 2 2$i2” from fitting the 
published KamLAND data 0. Only statistical errors are considered. Am^i has been marginalized in this plot. The 
vertical black line represents the best-fit sin 2 20 i 2 from the plane-wave analysis. 
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experimental support for such an assumption. In this paper we do not calculate or suggest the theoretical value of 
<j wp . We point out that the wave-packet impact is not significant for current reactor neutrino experiments. Nevertheless, 
the 1 a C.L. allowed range of cr wp is 0(0.1), within which the potential wave-packet impact could lead to significant 
effects and additional challenges in future neutrino oscillation experiments. 


3 Measuring neutrino mass hierarchy in reactor neutrino experiments 

3.1 Measurement of neutrino mass hierarchy 

The signs for Arris i and Arrisi have not yet been determined. Normal Hierarchy (NH) corresponds to positive 
Amsi and Am § 2 with v\ as the lightest mass state. Inverted Hierarchy (IH) corresponds to negative Arrisi and 
Arris 2 with V 3 as the lightest mass state [29 ] . 

As indicated by the recent data obtained by the Daya Bay Experiment, sin 2 20 i 3 = 0.084 ± 0.005 [26]. With such 
a (relatively) large value of $ 13 , it is possible to determine the neutrino mass hierarchy (MH) in reactor neutrino 
experiments at medium baseline (301 5 31 ] . 

For a detector at baseline L, the observed anti-electron neutrino flux of visible energy E v j s is given by 

f(E vis ) = I dEc/)(E)cr(E)Pee{L, E) 

■ R(E - 0.8MeV - E vis , 5E). ( 8 ) 

4>{E) is the reactor neutrino energy spectrum, and a(E ) is the inverse beta decay cross section. R is the detector 
response function with energy resolution 5E , which will be discussed in more detail in the next subsection. 

It is known that at a baseline of around 50 km, which corresponds to the first minimum of 612 oscillation for reactor 
neutrinos, the sensitivity for measuring MH is maximal. The upper panel in Fig. [7] shows f{E v i s ) at L = 53 km for 
NH and IH, with standard oscillation parameters in the plane-wave treatment (<r wp = 0 in Eq. ©)• 


3.2 Impact of wave-packet treatment 

Our wave-packet treatment shows that the amplitudes of neutrino oscillations will be reduced. In particular, damping 
of the 0 i 3 oscillations will be significant for an intermediate baseline reactor neutrino experiment if er wp is 0(0.1). The 
lower panel in Fig. [7] shows the neutrino visible energy spectrum at a baseline of 53 km, with standard neutrino mixing 
parameters and cr wp = 0.1. If <r wp is large, the neutrino spectra for NH and IH are indistinguishable from each other. 

We modify the GLoBES software |32l I33| to perform numerical simulations of a 53 km baseline reactor neutrino 
experiment, using a similar setup as in |34| and EQ: a 20 kton detector with 3% energy resolution, reactors with a 
total thermal power of 40 GW and a nominal running time of six years. In the absence of oscillations, the total number 
of events is about 10 6 . As this paper focuses on the wave-packet impact, the systematic errors were not taken into 
account in the following simulations 5 . We took the oscillation parameter values from global analysis [35] as Am 2 i = 
7.54 x 10 -5 eV 2 , (Ami 1 + Am§ 2 )/2 = 2.43 x 10 -3 eV 2 , sin 2 6 >i 2 = 0.307 and sin 2 0i 3 = 0.0241. 

To distinguish between NH and IH, we quantify the sensitivity of the MH measurement by employing the least- 
squares method, based on a y 2 function: 


,2 = (N z m - Nff 

^ / j jyM ’ 

i= 1 iV i 

where N ^ 1 is the measured neutrino events in the ith energy bin, and Nf is the predicted number of neutrino events 
with oscillations taken into account (without considering the systematic errors) 6 . The number of bins used -Wins is 
164, equally spaced between 1.8 and 10 MeV. 

We fit the hypothetical data set with Am 2 e as the free variable, defined as 

Z\m 2 e = cos 2 0 i 2 Am|! + sin 2 0 i 2 Am 2 2 . ( 10 ) 


The capability to resolve the mass hierarchy is then given by the difference between the minimum y 2 value for IH and 
NH: 


AyMH = niin(ym) - min(yN H )- 


( 11 ) 


5 We have also performed simulations with the following systematic errors m- 2% correlated reactor uncertainty, 0.8% uncorrelated 
reactor uncertainty, 1% spectral uncertainty, and 1% detector-related uncertainty. The results are similar to what we present here. 
6 Without loss of generality, in our simulations, NH is assumed to be the true mass hierarchy. The result is identical if we assume IH to 
be true. 
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Fig. 7: Upper panel: Expected neutrino visible energy spectrum for plane wave neutrino oscillation at 53 km for NH 
(solid) and IH (dashed). 

Lower panel: Same as in upper panel, but with er wp = 0.1. The two curves for NH and IH completely overlap due to 
the wave-packet impacts. 


Ax mh is used to explore the wave-packet effects as well as the impact of statistics and systematics in measuring 
the MH." If (j wp = 0, which corresponds to the plane-wave treatment, we get A\ mh ~ 19-5, implying that we could 
distinguish the MH with a confidence level of nearly 4 a, as shown in the upper panel of Fig. [8] However, if cr wp = 0.02, 
the sensitivity will be reduced due to the damping of oscillations and A\m h will drop to around 3.35, as shown in the 
bottom panel of Fig. [8] 

The solid (black) curve in Fig. [^further shows the variation of A\ mh as a function of <r wp . It shows that 
drops rapidly with cr wp , to become smaller than 3 a C.L. as er wp > 0.012. In this case, it is difficult to determine MH. 

In the rest of this section we will investigate how to increase the sensitivity to MH by improving the experimental 
setup. 


3.3 Consequences of energy resolution and statistics 


8 References [3611301137] suggest that the sensitivity of MH measurement depends on the detector resolution. As shown 
in Eq. (J8|, the observed z7 e flux depends on the detector response function R and energy resolution Se, which are 
defined as: 


R(E — Ey i S) 5E ) 



{E - E v i S ) 2 
2 (AE vis ) 2 1 


where the detector energy resolution is parameterized as 


5Ey i: 

E 


\ 


\[Ey is /MeV 


2 

+ h 2 . 


( 12 ) 


( 13 ) 


7 Note that min(x^ H ) and min(x 2 H ) can be located at different values of |Z\m 2 e |. 

8 The x-axes in most of the following figures just represent the U true-value” of Uwp* 



















9 
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Fig. 8: The sensitivities of measuring Arrive with and without wave-packet effects. NH is assumed to be the true 
hierarchy in the simulation. The solid (dashed) curve corresponds to the fitting of NH (IH, false hierarchy). o- wp is 
assumed to be 0 in the upper panel and 0.02 in the lower panel. 


In the previous subsection, we assumed that a = 0.03 and b = 0 in order to achieve a 3% detector resolution. 


As mentioned in the footnote in Section 2, the detector resolution is different from the energy uncertainty of 
detection at the microscopic level. The detector resolution is determined by the properties of the macroscopic detector, 
which should be taken into account even in the plane-wave assumption. Similar to the decoherence effect due to cr wp , 
a poor energy resolution can also destroy the measured oscillation effect. At L = 53 km, the $13 oscillation could be 
smeared out with a finite energy resolution, and the MH information could be destroyed, particularly in the case of 
large a wp - Fig. |10| shows that if <7 wp > 0.02, which is allowed by the Daya Bay and KamLAND data, the detector 


resolution (or the parameter a in Eq. 131, must be better than 3% in order to determine the MH with a C.L. of more 
than 2 a. 


On the other hand, we can also improve the MH sensitivity by collecting more data. Fig. |TT| shows the impact of 
statistics on the measurement of MH as a function of <7 wp , assuming a detector energy resolution of 3%, as suggested 
by references [5T1 i36j. Much longer run time would be required for a 2 a C.L. measurement if <r wp is larger than 0.01. 
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Fig. 9: Ax mh vs - ff w P - The experimental setup is described in the text. 



Fig. 10: Energy resolution required to resolve mass hierarchy at the 2 a (dotted), 3 a (dashed) and 4 a (solid) confidence 
levels, as a function of er W p. Other details of the experimental setup are described in the text. 


3.4 The optimal baseline 


Without considering the wave-packet impact, 50 ~ 60 km is the ideal location to measure mass hierarchy for reactor 
neutrino experiments [SI EE]. However, in the presence of significant wave-packet impact, longer baseline would lead 
to larger damping of the oscillation amplitude. Although reducing the baseline will lead to a loss of maximal phase 
difference between the NH and IH oscillation curves, this can save back part of the oscillation. Therefore, the optimal 
baseline of measuring Mil could be shorter than 50 km, depending on the value of er W p. 

shows the Axmh as a function of baseline for different values of <7 WP - In the case of the plane-wave limit 


Fig. 
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(d W p = 0), the MH can be distinguished with a confidence level of nearly 4 a at 53 km. However, as a, 
maximum Axm h shifts to shorter baseline. Fig. 


the 


13 


w P increases 

further shows the value of optimal baseline as a function of a wp - 
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cr W p 

Fig. 11: The number of neutrino events (no-oscillation) required to resolve MH at 2 a (dotted), 3 a (dashed) and 4 a 
(solid) C.L., as a function of <T wp . Other details of the experimental setup are described in the text. 


A x 1 



baseline (km) 


Fig. 12: A\ mh vs. baseline, assuming the true value of a wp = 0 (green dotted), 0.01 (red dashed), 0.015 (blue dot- 
dashed), 0.02 (black solid thick) and 0.025 (brown dotted thick) respectively. 


4 Conclusion 


The wave-packet impacts on current reactor neutrino oscillation and neutrino mass hierarchy measurements have 
been discussed. Our analyses show that the wave-packet treatment would not lead to significant modifications of the 
oscillation parameters from the Daya Bay and KamLAND results based on plane-wave assumptions. Moreover, our 
analyses also suggest that the energy uncertainty parameter a wp can be 0(0-1). 
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Fig. 13: The optimal baseline of measuring MH as a function of cr W p- 


The decoherence and dispersion effects depend not only on the initial neutrino energy uncertainty, but also the 
values of Am 2 and baseline. Since the measurement of the neutrino MH in medium baseline reactor neutrino exper¬ 
iments relies on the fast Am 2 2 oscillations, the decoherence and dispersion effects could be significant and make it 
more difficult. We found that even if <r wp is just around 0.02, the sensitivity of MH measurement would be largely 
reduced. The optimal baseline shifts to smaller value as a W p increases, due to the damping of oscillation amplitudes 
by wave-packet effects. 

We have to emphasize that we are not suggesting that the wave-packet impact would be so large as to make 
the mass hierarchy measurement in medium baseline reactor neutrino experiments impossible. We point out that the 
plane-wave model of neutrno oscillation is only an approximation, and the wave-packet treatment is more general. 
While the wave-packet effect is insignificant in the current reactor neutrino oscillation experiments, its impact on 
future oscillation experiments needs to be determined. 
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Appendix A: Derivation of the transition probabilities 

In order to derive the neutrino oscillation probability, we substitute Eq. Q into Eq. to obtain 
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2 i _, ■* 
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exp 


i(p v z — Ei(p v )t) - 
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(u\vp F/(/+)) (z Vi(pis)t)^ 

2 7TTT 
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(A.l) 


■i), (A.2) 


G v Gv 

where a wp = —r « -= 7 —r, 

Ei(p v ) E(p u ) 


C Pl+m 2 f /2 _ Ei(p u ) 3 ~ E{p„) 


which is a moving Gaussian wave packet with dispersion. Here, we have used the approximation Ei(p u ) ss p v ss E(p„), 
where E(p u ) is the neutrino energy measured by the detector. Eq. (A.2) shows that the propagating state | i/i(z,t)) 
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exhibits dispersion - a spreading of the wave packet in space. Both the group velocity and the rate of dispersion depend 
on the neutrino mass. The heavier the mass, the smaller is the group velocity and the higher is the rate of dispersion. 

Let \v a ) (a = e, (i,.. .) represent the neutrino flavor states. Since the flavor eigenstates are superpositions of mass 
eigenstates \vi), from Eq. (|A.2[), the time evolution of a flavor state is given by 


\v a (z,t)) = 'Y^Uai\vi{z,t)) 


=£ 


u* 

rv 


l + i- 


?E(p„) 


exp i(p v z - Ei(p v )t) - 


(^wpl?(pi/)) (z Vi (pu)t) 


I £ UpiWp) ■ 

\P=e,iJ,... J 

Then the transition probability of v a —> vp at a distance L and time t is given by 


(A.3) 
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f 2 a \\ 
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V T i T J J J 


(A.4) 


which is a function of both time (t) and distance (L). 

In an oscillation experiment, the neutrino is detected at a fixed baseline L but the time is not measured. In order 


to obtain the oscillation probability as a function of the baseline, the time has to be integrated out in Eq. (A.4). Since 


reactor neutrinos propagate almost at the speed of light, P a p(L,t) is non-zero only around t ~ L. The transition 
probability is non-zero only within a small time window At around a time where 


L= v 1 + v^ 
2 


— t £z, — 


2 L 

Vi + Vj 


L. 


(A.5) 


On the other hand, the size of At is constrained by the spatial width of the wave packet which is typically much 
smaller than the baseline L , which means that At -C tz,. Moreover, 


3(1 


1 

Ti 


dt E(p„) " wp 

d(z - Vi(jp v )t) 


dt 


= \Vi(pu)\- 


(A.6) 

(A.7) 


Eq. (A.6) shows that the factor (£/ty) changes slowly with the variable t. Within a small At, these terms can be treated 
as constants. Therefore, within the small integration region which is constrained by the width of the wave packet, it 
is acceptable to approximate t = th ~ 2 for this factor, 


2 2 
^ w 

E{p v 


t rrii cr wp t ?n^(T W p 


E(p v ) 


(A.8) 


However, since Eq. (A.7) shows that the derivative of (z — Vi(p^)t) is not negligible even in a small region. We did not 
use the same approximation in the factors (z — Vi{p v )t ) and (z — Vj{p u )t). Therefore, the integral of Eq. (A.4) can be 
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approximated as, 


■xp{z) = I P a p(z,t)dt 
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In the last step of Eq. (A.9), we have ignored the terms proportional to 

nej 
by 


EM 2 

negligible for ultrarelativistic neutrinos. In the case of reactor neutrino experiments, the v e survival probability is given 


m 2 m 2 


(A.9) 


, since they are expected to be 
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(A.10) 


The damping factor exp(—7^ ) corresponds to an decoherence effect from delocalization, which is neglected in the main 
article. This is because in most circumstances, 7^ ~ 0. The details of 77 and the delocalization decoherence effect 
will be discussed in Appendix B. 


Appendix B: The decoherence effect from delocalization 

The decoherence effect discussed in Section 2 is due to the separation of different neutrino wave packets. With larger 
values of er wp , the decoherence effect would be more significant. On the other hand, there exists another decoherence 
effect which is due to the delocalization of the production and detection processes. Different from what we have studied 
above, the decoherence effect from delocalization will become significant only if <7 wp is extremely small, or the spatial 
uncertainty cr x is large. In fact, in neutrino oscillations, one of the coherence conditions is that the intrinsic production 
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(and also detection) energy uncertainties are much larger than the energy differences between different mass eigenstates 
AEij = Ei — Ej (Ei is the energy of mass eigenstate \vi)) [38], namely, 


AEu = Ei - Ej 


Arrij 

E v 


o'i* — EvOv 


which means o x -C L°| c (ov ~ 1 /o x ). 


(B.ll) 


Eq. ( B. 1 ip implies that in order to measure the interferences between different mass eigenstates, the spatial uncertainty 
o x , has to be much smaller than the oscillation length L°J C . 

Eq. 0 is just an approximate neutrino oscillation probability formula and it does not describe the decoherence 
effect from delocalization. More precisely, the flavor transition probability in Eq. © should be multiplied by an 
additional factor 


ex P(—7ij), 


(B.12) 


where 


7 ij — 


(L 


°!j c ) 2 E 2 c 


(i + yfi) (£° sc ) 2 ’ 


With these delocalization terms also taken into account, a more complete v e survival probability should be written 


as Eq. (A. 10> in the previous section. The damping factor in Eq. (B.12) corresponds to the decoherence effect from 
0, the modifications from delocalization are negligible. In this case Eq. (A.101 just reduces to 

rr 2 - 

T osc 
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delocalization. If 
Ecp (5). 7 ij is proportional to 


( T ° sc ) 2 


and so the decoherence effect from delocalization matters only when o x 


about 1 km in the measurement of 0\z oscillation. In this case o wp is extremely small. 

If we use Eq. (A. 101 rather than Eq. © to do analysis, we will find that the 7 ij term will offer a lower bound on 
CT W p and the delocalization effect is significant only when o wp is extremely small, which implies a large o x . Fig. [14] 
shows our analysis of the Daya Bay data, which is similar to Fig. [5] but this time the delocalization term is considered. 


Marginalized over A m 31 2 



sin 2 (2f? 13 ) 

Fig. 14: Same as Fig. pi but delocalization is taken into account and the y-axis is in log scale. 


From Fig. [14] we can see that in a large range of the parameter space, there are no modifications on the true value 
of sin 2 20i3. This means that in most of the parameter region, wave-packet impact can be safely neglected. Moreover, 
Fig. 14 also suggests that only if o wp ~ O(10 -16 ), which means that o x ~ 0(1 km), the decoherence effect from 


delocalization is significant. 

At this point, we can discuss the wave-packet impact in two different regimes. If <r W p is large ( 0(10 x )), since Ai 


oc <7„ 


the damping factors exp(— \ij) in Eq. (1 51 become significant and the decoherence effect from the separation of 


wave packets cannot be neglected. On the other hand, if < 


0(10 °) or even smaller, then the additional damping 


factor in Eq. (B.12 I starts to dominate since 7 ij oc o wp . In this case the decoherence effect from delocalization becomes 
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important. Nevertheless, in most reactor neutrino experiments, the dimensions of the reactor cores and detectors are 
just around a few meters. It is unlikely that the spatial width of the initial neutrino wave packet would be larger than 
1 km. 

However, in most reactor experiments, including current ones such as Daya Bay and KamLAND, and also the 
proposed measurements of neutrino mass hierarchy at medium baseline, the delocalization terms exp(—7,,-) can be 


neglected. This is because 


<C 1 (or, the spatial uncertainty a x <C L°| c ). In this case, Eq. 


(L°f) 2 £% 2 p 

reduces to Eq. ©. Therefore, we neglected the decoherence effect from delocalization in our study and performed the 
simulations and analyses with Eq. ©■ 


B.12) just 
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